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Abstract. We consider a class nonlinear Schrodinger / Gross-Pitacvskii equations (NLS/GP) with a 
focusing (attractive) nonlinear potential and symmetric double well linear potential. NLS/GP plays a 
central role in the modeling of nonlinear optical and mean-field quantum many-body phenomena. It is 
known that there is a critical L 2 norm (optical power / particle number) at which there is a symmetry 
breaking bifurcation of the ground state. We study the rich dynamical behavior near the symmetry 
breaking point. The source of this behavior in the full Hamiltonian PDE is related to the dynamics of a 
finite-dimensional Hamiltonian reduction. We derive this reduction, analyze a part of its phase space and 
prove a shadowing theorem on the persistence of solutions, with oscillating mass-transport between wells, 
on very long, but finite, time scales within the full NLS/GP. The infinite time dynamics for NLS/GP 
are expected to depart, from the finite dimensional reduction, due to resonant coupling of discrete and 
continuum / radiation modes. 



1. Introduction and Outline 

The cubic nonlinear Schrodinger / Gross-Pitacvskii (NLS / GP) equation 
(1.1) id t u=(-A + V(x))u -\u\ 2 u 

plays a central role in the mathematical description of nonlinear optical and quantum many-body phenom- 
ena. In the context of nonlinear optics, u denotes the slowly varying envelope of a nearly monochromatic 
electromagnetic field propagating in a wave-guide, t the distance along the wave-guide and x £ K 2 the 
dimensions transverse to the wave guide [5] ■ At low intensity, light is guided to a higher refractive 
index region, corresponding to a potential well V(x). The Kerr nonlinear effect gives rise to an increase 
in the refractive index in regions of higher intensity, \u\ 2 , and therefore a "deeper" effective potential well 
V(x) — \u\ 2 . In the context of quantum many-body physics, NLS/GP emerges in the mean-field limit of 
many weakly interacting identical quantum particles obeying Bose statistics, as the number of particles 
tends to infinity pM [TtJl |3"0] and other recent works. The potential V(x) governs the confining trap for the 
quantum particles □ 

In this paper we focus on a class of symmetric double-well potentials. A model to keep in mind is the 
two parameter family of symmetric double well potentials: 



(1.3) V(x) = V(x; a, L) = - 



1 (m - L ) 2 1 (g±jOj 

& An-* -J- & drr'^ 



V 47TCT 2 V47T<7^ 



r 2 

which converges as a — > + to a double-delta well at ±L. In figure [T] the case L = 3, a = 1 is shown. 



^NLS / GP falls into a larger class of models: 
(1.2) id t u = (-A + V(x))u + gK^\u\ 2 Ju, 

allowing for more general nonlinear terms, e.g. more general functions of \u\ 2 or K nonlocal; see, for example, |21| . For g < 
(focusing case) our analysis yields very similar results to those obtained above. For g > (defocusing) our methods can be 
adapted to prove a shadowing result for trajectories of the finite dimensional reduction which arises. 
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Figure 1. Double-well potential, V(x; 3, 1). Superposed are the two eigenstates, even 
ground (^o) and odd excited {tpi) states, of the Schrodinger operator, —d 2 + V(x; 3, 1) 



Double-well potentials are of particular interest in optics as models of coupled parallel wave guides 
channeling light, which interacts through through evanescent tails. In the quantum context, these are 
natural and simple systems in which to study quantum tunneling. As discussed below, the combined 
effects of a confining double-well potential with focusing cubic nonlincarity lead to the phenomenon of 
spontaneous symmetry breaking of the ground state at sufficiently high optical power or particle number. 
See, for example, |4Q[ [4j I19j for experimental studies of symmetry breaking. Our goal in this paper is to 
investigate the phase space dynamics of NLS/ GP near the symmetry breaking point. 

The equation NLS/GP, (|1.2| . is a Hamiltonian system and expressible in the form: 



(1.4) 



■a 5H 



V\u\ 2 --\u\ 4 



dx. 



Ti. denotes the conserved Hamiltonian energy functional: 

H[u] = J 

The conserved squared L 2 norm (particle number / optical power) is denoted: 
(1.5) N[u] 



\u\ 2 dx. 



We are interested in the dynamics near special classes of nonlinear bound states of NLS / GP. Nonlinear 
bound states arc solutions of the form 

u(x,t) = e- mt tfnOr), 

where is spatially localized: 

(1.6) (-A + V(s))*n - |*n| 2 *n = fi *n, *neH\R). 
Consider first the linear case of the linear eigenvalue problem: 

(1.7) (-A + V(x))*n = VneH^R). 

In this case, there is a least energy ground state, ipo with corresponding simple eigenvalue f^o [29]. If 
the separation between wells is sufficiently large, then the ground state eigenfunction is a bimodal positive 
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Figure 2. A symmetric (bimodal) soliton and an asymmetric soliton, which exists only 
for Af > N cr - 

symmetric state, which reflects the discrete symmetry of the potential [El EH E] ; see figure [TJ In addition, 
there is an anti-symmetric (odd) state, tpi with energy f2i, such that ilo < fii < 0. 

In the attractive / focusing nonlinear case, the character of solutions, and the solution set varies 

with the solution norm. Indeed, if we consider the set of solutions of (|1.6p on the level set 

(1-8) J\* n \ 2 =Af, 

we find that for large enough well-separation, there is a symmetry breaking threshold Af cr [21| ; see figure^ 

(1) If Af < Af cr there is a unique positive, symmetric and bimodal state. 

(2) For Af > Af cr , (modulo phase) there are three positive localized states: a symmetric state (which 
exists for all Af > 0) and two are asymmetric states, biased respectively to the right and left wells. 

(3) As Af increases beyond Af cr this symmetry broken state becomes increasingly concentrated in one 
of the wells [11 [21]. 

(4) The symmetric (bimodal) state is dynamically stable for Af < Af cr and unstable for Af > Af C r- For 
Af > Af cr the asymmetric states are stable. 

That symmetry breaking occurs at sufficiently large values of Af was studied variationally in [7] for the 
nonlinear Hartree equation. Their method can be adapted to a large class of equations, for which a ground 
state can be realized as a minimizer of a Hamiltonian, Ti., subject to fixed Af. The above detailed portrait 
of the symmetry breaking transition and the exchange of stability among branches was studied in detail in 
|21j . In particular, for large well-separation, L, 

(1.9) AA cr (L) = O(fi 1 (L)-fi (i))=C'(e- KL ), k > 0. 

Figure [3] displays a numerical computation of the symmetry- broken state, occurring for Af > Af cr , where 
V(x) is of the type displayed (|1.3[) . 

Our goal is to explore the detailed general dynamics near the symmetry breaking transition. Toward 
a formulation of precise results, wc first introduce a class of double-well potentials in dimension one 0. 
Following [12], start with a single rapidly decaying potential well centered at 0, Vo(x), for which the 
Schrodiner operator H = — d% + Vq has exactly one (simple) eigenvalue u. Then, construct a double well 



The results of this paper can be extended to higher dimensions. 
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Figure 3. A numerical plot of the symmetry breaking in the soliton curve. Specifically, 
we plot here N versus fl for a double- (Gaussian) well potential with separation parameter 
L = 3 and a = 1. Note the stable curve is invariant under reflection about the y axis, 
hence the slightly thicker curve represents the resulting degeneracy for the asymmetric 
states. 



There exists L > such that for L > L , has a pair of simple eigenvalues, f2 = ^ (L) and = Sl\{L) 
and corresponding eigcnfunctions (even) and ■01 (odd): 



As noted above, the symmetry breaking threshold, N cr (L) (equation (|1.9p ) is exponentially small for large 
well-separation. Therefore, to study the dynamics in a neighborhood of the symmetry breaking point, it 
is natural use coordinates associated with the linear operator Hq = —d^ + V. Throughout this result, 
we will assume that V is such that ipo and are the only discrete eigenfunctions of Hq and in addition 
that V is sufficiently smooth and decaying as defined in [39] and discussed in Appendix [B] as to guarantee 
dispersive estimates to do perturbation theory. We note that these assumptions will be satisfied by double 
delta wells, for which we have done our numerical computations. These ideas will be explored further in a 
forthcoming note [S] . 

Thus, we expand the solution as follows: 




V {x- L) + V {x + L), L>0 



Q < ^i < 0. 



(1.12) 



u(x,t) = co(i)^o(a;) + ci(t)ip 1 (x) +R(t,x), 
(^,R(-,t) ) = 0, j = 0,l. 
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By orthogonality, we have 

(1.13) MM;t)] = M*)| 2 + |ci(t)| 2 + f \R(x,t)\ 2 dx = A/>(-,0)]. 

Referring to this decomposition, we now give an overview of the paper: 

(1) In section [2] we express the NLS / GP as an equivalent dynamical system governing co(t), C\(t) 
and R{x,t). This Hamiltonian system, equivalent to NLS / GP, has the form of two equations 
governing discrete nonlinear "oscillators", coupled to an equation for a wave field, R(x,t). 

(2) In section [3] we study the finite dimensional reduction governing cq and c\, obtained by dropping 
all terms which couple the oscillator and field variables. This reduction is a finite dimensional 
Hamiltonian system with conserved Hamiltonian: 

-ff(c ,c ,ci,ci), 

given in (|3.3p . and I 2 invariant 

Aqc ,ci] = |c | 2 + | Cl | 2 . 

Remark 1.1. Note also that the reduction above is what we get if we make the Ansatz (|1.12|) into 
the Lagrangian of NLS/GP, that is restrict Cnls-gp to the bound state manifold, and obtain the 
equations of motion as in [15] . 

Use of these two invariants facilitates an analysis of the finite dimensional phase space; trajec- 
tories with prescribed values of H and N lie on a two-dimensional surface. We then discuss some 
of the rich dynamics of this finite dimensional reduction and our goal is to prove their persistence, 
for non-trivial time scales , for with the full PDE, NLS/GP. 

For the dynamics on level set TV ~ N^ r D as N cr , we establish the following behavior: 

(a) N < N^ r D : There is an elliptic fixed point, corresponding to the stable symmetric state of 
NLS/GP for N < N cr . A ncighorhood of this fixed point is foliated by stable time-periodic 
solutions. 

(b) N > N£P : The fixed point for N < N^ r D persists, but transitions from being a stable elliptic 
point to an unstable saddle. This corresponds to the unstable symmetric-bimodal state for 
Af > M cr - At N = N^ r D . two new elliptic equilibria bifurcate and, a neighborhood of each is 
foliated by stable time-periodic solutions. These stable time-periodic oscillations correspond 
to stable oscillations around the stable asymmetric standing waves for M > JV cr . 

(c) N > N^ D : There are periodic solutions, outside a separatrix, which encircle both new 
equilibria. In the physical configuration space, these correspond to soliton transport from 
one well to the other and back, continuing periodically. Furthermore, one can quantify the 
"energy barrier" that must be exceeded to dislodge a "soliton" from localization about one of 
the wells. Energy thresholds, such as that described above play an important role in transport 
of energy in inhomogeneous and discrete systems. Natural directions to pursue beyond this 
work are transport in systems with many wells and the Peireles-Nabbaro barrier for motion of 
localized coherent structures discrete lattice systems. 

(3) In section [5] we prove that certain phase space structures for the finite dimensional dynamical 
systems, persist in character for the full NLS/GP system, on long, but finite time scales. Stated 
nontcchnically, 
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Theorem: For any sufficiently small amplitude periodic solution about equilibrium an equilibrium 
state of the finite dimensional reduction (N above or below the bifurcation threshold), there is a 
solution of the PDE (NLS/GP), whose projection into the finite dimensional phase space, shadows 
this finite dimension orbit on very long time scales. 

A precise statement is given in Theorem [T] sec also Figure 01 The time scales on which these 
results hold enable us to see nearly periodic oscillations for the PDE on long time scales, i.e. 
through many, many oscillations, but not on an infinite time scale. Indeed, we do not expect the 
persistence of such oscillations on infinite time scales due to nonlinear coupling of bound to to 
radiation modes for the full system, for example, [351 HI] ■ We conjecture that in the infinite time 
limit, the soliton executes a (radiation) damped oscillation to some stable nonlinear bound state, 
corresponding to the damped oscillatory decay to a stable equilibrium of the finite dimensional 
reduction. Evidence is presented in section [6l where numerical simulations are discussed. 

We remark that our current theorem does not apply to perturbations of general "large" periodic 
orbits of the finite dimensional reduction. More detailed information on the Floquet theory of the 
linearized equations about general periodic orbits is still needed. This is currently being investi- 
gated. 

(4) In section [7] wc provide a summary and discussion of open problems. 

Acknowledgments. JLM was supported, in part, by a U.S. National Science Foundation Postdoctoral 
Fellowship and a Hausdorff Center Postdoctoral Fellowship. MIW was supported, in part, by U.S. NSF 
Grants DMS-04-12305 and DMS-07-07850. MIW wishes to thank Vered Rom-Kedar and Eli Shlizerman 
for stimulating discussions. 

1.1. Notation. 

(1) The spaces iP(R n ), L p (R n ), W k ' p are the standardly defined Sobolev integration spaces. 

(2) We have the L 2 inner product: (/, g) = J R „ fg. 

(3) Projection onto the bound states of —A + V: 



2. Formulation of NLS/GP as a coupled finite-infinite dimensional system 

In this section we derive an equivalent formulation of NLS / GP, appropriate for studying the exchange 
of energy between the bound and radiative parts of the solution. We substitute the decomposition (|1.12p 
into NLS / GP and, to the resulting equation, apply the projection operators Po, Pi and P c to obtain the 
coupled system 




(irjf)^, j = 0,1. 



(2.1) 



ic - fi c + a 00 oo|co| 2 co + aooii(c 2 c + 2\c 1 \ 2 c Q ) = F (co, Ci, Co, ci; R, R), 
ici - QxCi + ami|ci| 2 Ci + aoon(coCi + 2|c | 2 ci) = F 1 (c ,ci,c ,_ci;R,R), 
iR t - HR + P c F b (c Q ,ci,c ,ci) = P c F R (c ,ci,c ,ci;R,R). 



Here, 
(2.2) 



<Hjkl = {ipi^jtpkjtpi) 
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Figure 4. A numerical plot of an oscillatory solution to at times a, b, c, and d 

respectively from the corresponding phase plane diagram of the finite dimensional Hamil- 
tonian truncation. 



(2.3) Fj = tt 3 F, 
for j = 0,1 where 

F = [2|c | 2 ^o+2|ci |Vi + 2(coCi + ciCoty'oV'i] # 

(2.4) + [ c^i + CoV^-o] # 2 + [2c ^o + 2 Cl V^i ] |i?| 2 + li?! 2 ^, 
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and 

(2.5) F b = P c [ |c | 2 c ^ + (egg! + 2|c | 2 c )^V'i + (c 2 c + 2c \ Cl \ 2 )Ml + N'W? ] , 

(2.6) F R = P c ( [ 2|c | 2 V'o + 2|ci| 2 V> 2 + 2(c ci + CiCo^l ] R 

+ [ c 2 ^l + c\ip\ + 2c cit/;oV'i ] R 

+ [ c^i + c ?Ao] R 2 + [2coVo + 2ciV>i ] |i?| 2 

+ \R\ 2 R). 

We have that a^fc; ^ ^=^> i + j + k + I = mod2. For simplicity we take 

a.ijki = lifi+j + fc + / = mod 2. 

Initial conditions for (|2.1[) are: 

(2.7) Cj (0) = (^,<-,0)), j = 0,l; i?(-,0)=P c i?(-,0). 

The system (|2.ip may be viewed as an infinite dimensional Hamiltonian system, comprised of two coupled 
subsystems: one finite dimensional, governing bound state degrees of freedom described by Co and c\, and 
a second, infinite dimensional, governing a dispersive wave field, i?(-, t) = P c i?(-, t). In the next section, we 
focus on the finite dimensional truncation of (|2.ip obtained by setting R — 0. After obtaining a detailed 
description of the phase space of this truncation in section [3j we then turn toward proving the long-time 
persistence of structures within the finite dimensional system, within the full infinite dimensional problem. 

Before embarking on this path, we conclude this section with an alternative coordinate description of 
(|2.ip . These coordinates prove very useful in understanding the bifurcations within the finite dynamical 
subsystem, as Af is varied, and its participation within the infinite dimensional dynamics. 

2.1. Alternative coordinates. We shall require the following decomposition, proved in Appendix [A"l 
Introduce the following change of coordinates: 

( c (0,ci(t),c7,(i),cl(i),ii(-,t) ) i * ( A(t),a(t),(3(t),6(t),R(-,t) ) 

defined by 

(2.8) o(t) = A(t)e^\ 

(2.9) Cl {t) = {a{t)+if5{t)y e{t \ 

(2.10) u(x, t) = e w(t) (A(t)V>o + («(*) + iP(t))ipi + R). 

Such coordinates have been used for finite dimensional systems in for instance the works |23j , |33j . 
Substitution of this ansatz into NLS/GP, we see 

(iA — OA — Sl A)il) + (ia-P- 6{a + i/3) - (a + i/3)fii)^i 
+iRt - HR-6R = F(A, a, (3; R, R), 

where F is determined as in (|2.4[) . Hence, 

A = lm(w F), 

a = ^/3 + Qi/3 + Im(7ri(F)), 

$ = -6>a-fiia-Rc(7Ti(F)), 

= -no-^Re^o^)), 
iR t = HR + 9R + P C (F) 
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or 

A = Im(7r F), 

a = (-n -A- 1 n (F))(3 + n 1 (3 + lm(TT 1 (F)), 

[3 = -(-fi -yl" 1 Re(7ro(F)))a-fi 1 a-Re(7ri(F)), 

= -f2 -^" 1 Re(7r ( J F)), 

iRt = (H -n )R-A- 1 Rc{MF))R + Pc(F). 

Note, there will be linear terms in R and R contained in F, which we must handle very carefully in our 
eventual iteration argument. 

This leads to the system of the form 

A = -2apA + Errov{R,R;A,a,/3), 

a = [Oi - (a 2 +/3 2 ) -A 2 + 9](3 + Ettot(R, R; A,a, f3), 

$ = -[tt-L - (a 2 +0 1 +A 2 ) - 2A 2 + % + Error A, a, 13), 

6 = -n + A 2 + (3a 2 +f3 2 ) + A^ErmriR, R; A, a, (3), 

iR t = (H — fl )R — A~ x Re(ir (F))R + P c F b (A, a, (3) + P C F R (A, a, (3; R, R). 

Then, substituting 9, the system becomes 

(2.11) A = -2a(3A + Error A (R,R;A,a,/3), 

(2.12) a = -n + 2a 2 }(3 + ETTor a (R,R;A,a,(3), 

(2.13) $ = -[Q 1 -n a -2A 2 + 2a 2 ]a + EiTOi /3 (R,R;A,a,P), 

(2.14) 9 = -too + A 2 + 3a 2 + (3 2 +Ettot 9 (R, R; A, a, (3), 

(2.15) iR t = (H -n Q )R+{A 2 + 3a 2 + (3 2 )R + P c F b (A, a, (3) + P C F R (A, a, (3; R, R), 

with P C R = R, 

For simplicity, set a = (a, (3, A). 

Proposition 2.1. Assume A, . . . ,R defined as in (|2. 1 1[> - (|2. 15() . Then, we have 



(2.16) 


| Error a \ 


< 


I^ a lli2iu» + i^iiiiiii» + iiiiiii», 


(2.17) 


| Error a \ 


< 


\<t\ 2 \\R\\ l ~ + \<t\\\R\\U + \\R\\U, 


(2.18) 


| Error p \ 


< 


|ff| 3 Hii|iL- + |ff|iii2|ii- + iiij|ii- ) 


(2.19) 


| Errors \ 


< 


A" 1 (\ff\ 2 \\R\\ L ^ + 1^1^ + \\R\\l~), 


(2.20) 


\P c F h \ 




\P C [A 3 ^ + (a 2 + f3 2 ){a + iflift + A{a + if3) 2 ^ 






+ 


2A(a 2 + /3 2 )V>jVo + A 2 (a - ^/?)^/> 2 V'l + 2A 2 {a + i0)^i] \ 


(2.21) 


\PcF R \ 




\P C ([2A 2 i> 2 Q + 4A(a)^i + 2(a 2 + (3 2 )^ 2 } R 



+ [AVo + {a + if3) 2 i>l + 2A{a + i(3)^i] R 

+ [Ayj Q + (a - 0)^} R 2 + [2Ayj Q + 2{a + if3)^] \R\ 2 + \R\ 2 R) \. 

The error terms in a and (3 gained from 9 are of the form A~ 1 a, A~ x (3 <C 1. 

Proof. These follow directly from the computations in Appendix [A] □ 
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(2.22) 



For convenience, we rewrite the system (|2.11[) - (|2.15[) in compact form as: 

cr = F FD (a)+G FD (a;R,R), 
= -fi + A 2 + 3a 2 + 1 + G 9 (R, R; A, a, (3), 
(id t - (H - Q Q ) - (A 2 + 3a 2 + I3 2 ))R = P c F b (a) + P c F R (a; R, R). 

For simplicity, we take F b = P c F b and Fr = P c Fr in the sequel. 

3. Phase space of the finite dimensional Hamiltonian truncation of NLS/GP 

In this section we study the finite-dimensional system obtained by setting the dispersive part of the 
solution, R(x, t), equal to zero in (|2.1|) . We denote the solution of the resulting system by (po(t), pi(t)) £ C 2 : 

ipo = Qopo - (popo + 1p\p\pQ + plpo), 
ipo = -OoPo + (p~o 2 Po + ZpiP'iPo + Px 2 Po), 

ip\ = Qi/9i - (pfpt + 2p poPi + PoPx), 
ipi = -Clipi + (pi 2 pi + 2\po\ 2 p~i + Po 2 Pi)- 

Symmetry breaking for a system of this type arising from a general class of defocusing nonlincaritics was 
considered recently in [31] . 

This is a two degree of freedom Hamiltonian system, with time-translation and phase invariances inher- 
ited from NLS / GP. The associated time-conserved Hamiltonian and L 2 (optical power or particle number) 
functionals are: 

(3.2) H = no\ P o\ 2 + fiiM 2 - \\pof - ^|pi| 4 ~ 2\ Pl \ 2 \p \ 2 - \{pW + pi 2 pg), 

(3.3) N=\ Po \ 2 + \ Pl \ 2 . 

In terms of H , system (|3.1|) can be expressed in Hamiltonian form 



(3-1) 



(3.4) 
where 

(3.5) 



idtp = JVpH, 



Po 




Po 


, J = 


Pi 




. Pi . 





10 
-10 








1 

-1 



ft. 



In terms of the alternative coordinates of section (12. ip , 

p a = Ae l9 , pi = (a + ij3)e v 

hence the system (|3.ip takes the form: 

a = \ fii +0- {a 2 + (3 2 ) -A 2 

) $ = -[si 1 +0-(a 2 +/3 2 )-3A 2 j a, 
A = -2aPA, 
= -n o + A 2 +3a 2 +P 2 . 

Recall that, for simplicity, we have set OijU = 1 for all i, j, k, I = 0, 1. Note that 9 completely decouples 
from the A, a, (3 equations, meaning wc have 

a=[n 1 -n Q + 2a 2 ] /3, 
(3.7) { = - [ fti - n Q - 2A 2 + 2a 2 ] a, 

A = -2af3A, 
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where now the 9 equation is decoupled to give 

(3.8) 9 = -O + A 2 + 3a 2 + 2 . 

One can verify changing coordinates in (|3.2II3.3|) . or directly from (|3.6[) . the time-conserved quantities: 

(3.9) N = A 2 +a 2 +p 2 , 

(3.10) H = n A 2 + n 1 {a 2 + p 2 ) - ^A 4 -h,a 2 + p 2 ) 2 -2A 2 (a 2 + f3 2 ) - A 2 (a 2 - p 2 ). 

We obtain a closed system for (a, (3) using that N = A 2 + a 2 + p 2 is conserved. Then, we may reduce 
the system to 

(3.11) d= [O 1 -fio + 2a 2 ] 0, 

= - [ Oi - O + 2(2a 2 + P 2 - N) ] a. 
In addition, the system has the conserved quantity 

H = N{^-^) + ^ w {a 2 + f3 2 )-2Na 2 + {a 2 +l3 2 )+a i -l3 i . 

Now, let us define the matrices B eq and B eq to be those related to linearization about the equilib- 
rium solution (A eq , a eq , /3 eq , 9 eq ) and the decoupled reduced system {A eq , a eq , /3 eq ) respectively. Simi- 
larly, let M(t) and M be the resulting monodromy matrices for nearby time dependent periodic orbits 
(A(t),a(t), P(t),9(t)) and the decoupled reduced system (A(t),a(t), f3(t)) respectively. 

In the following section, we actually discuss the relevant bounds on the operators e B "' t and e Bt . 

3.1. Bifurcation of equilibria for (|3.7p and Symmetry Breaking in NLS/GP. Recall that N = 
N[A, a, (3] is a constant of the motion for the system (|3.6j) . corresponding to the physical quantities optical 
power or particle number. Thus it is natural to explore the nature of the phase space restricted to the level 
sets of N. We are interested in time-periodic states of frequency fl, corresponding in the physical space to 
solutions of NLS/GP of the form u(x, t) = e~ lflt U. Thus, we transform the system to a rotating frame by 

9(t) = 9(f) - fit 

fix - ft + Q(t) - (a 2 + /3 2 ) - A 2 J /?, 

O x - + 0(t) - (a 2 + 1 ) - ZA 2 1 a, 
-2apA, 

l-n Q + A 2 +3a 2 +f3 2 . 
The states we seek are equilibria in this rotating frame. Thus, we have 

(3.17) [ Oi - O - (a 2 + 2 ) - A 2 ] = 0, 

(3.18) [ Oi - O - (a 2 +0 1 ) - 3A 2 ] a = 0, 

(3.19) af3A = 0, 

(3.20) - O + A 2 + 3a 2 + [3 2 = 0, 
whose solutions we consider on the level set 

(3.21) A 2 + a 2 + p 2 = N. 

11 



setting 




(3.12) 




and obtain 




(3.13) 


Q 


(3.14) 





(3.15) 


A 


(3.16) 


e 
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It is easy to observe an equilibrium corresponding to the 
Symmetric states: 

(3.22) For7V>0: A, = N? , a* = /3* = 0, = fl - N. 
Via (|2.10[) wc identify this with the symmetric ground state of NLS/GP; 

u{x,t) ~ e *(-n«+N)t 

Due to its correspondence with the symmetric state, we refer to this equilibrium of the finite dimensional 
reduction, the symmetric equilibrium. 

A second, bifurcating family can be found explicitly as follows. Define 

(3.23) NZ? = ^ = °LZ°°. 

Remark 3.1. Had we set the nonlinearity coefficient g = —1 and the "interaction weights", a^u equal to 
one, we would have 

(3 24) N FD = — > 0- 

which is easily observed by comparison to a single well potential symmetric state for L sufficiently large 
(see EH " 



Using (j3~2TJ)l we first eliminate from (|3~T7)) and obtain [ ^ - f2 + 2a 2 ] (3 = 0. 
Since fix — J7o > wc conclude (3 = 0. Thus. (|3 . 1 9[) is satisfied. We now use (|3.20[) again to eliminate f2 
from (|3. 18|) . Thus we have, since (3 = 



A 2 - a 2 = N™, 
A 2 + a 2 = N. 



Solving for A and a we have the following equilibria, corresponding to symmetry broken states, which 
bifurcate for at N = Nl. D : 



Symmetry broken states: 

For N > NZ D = 



rtio 



' N 4- N FD \ 1 fN-N FD \i 
(3.25) A* = ( — — — — — — — , a* = , & = 0, fi* = -2Q Q + - N. 



Via (|2.10j) wc identify this with the asymmetric ground state of NLS/GP; 

Due to its correspondence with the asymmetric (symmetry-broken) states of NLS/GP, we refer to this 
equilibria of the finite dimensional reduction as asymmetric equilibria. 



■^Though N F r D is a good first order approximation to N cr a, the true symmetry breaking point in the nonlinear problem 
11. H . As in this note we will be studying existence of solutions close to those described by the finite dimensional dynamcis 
of 13.11 1 , wc will work from here on using N cr = N F r D . 
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(3.26) 



3.2. Stability of equilibria; finite dimensional analysis. We consider the stability of the various 
solution branches obtained in the previous section. We rewrite the system (|3.13|) . using the last equation 
to eliminate —51 + from the equations for a and (3. Thus we have 

a = [ Oio + 2a 2 ] ft 
/? = - [ fiio + 2a 2 - 2A 2 ] a, 

A = -2a(3A, 
<d = n-fl a + A 2 + 3a 2 + (3 2 . 

Note that in these coordinates the equations for a, (3 and A decouple from the equation for 0. 

The finite dimensional system has a phase portrait, equivalent to (|3.1|) . sec Figure [5] In particular, we 
observe elliptic and hyperbolic equilibria, and periodic orbits. We now embark on detailed linear stability 
analysis of these states. 

Linearization about an arbitrary solution 

(3.27) (a*(t),/3*(t),A*(t)A(*)) 
gives the linearized perturbation equation 



(3.28) 





5a 




4a*ft 


(fiio + 2a 2 ) 





" 




8a 




5(3 




-(Oi + 6a 2 -2A 2 ) 





-Aa,A, 







8(3 


SA 




-2 A, p, 


-2a*A* 


2a,P, 







8A 




59 




6a* 


213, 


2A, 







88 



(3.29) 



B(t) 



5 a 
5(3 
8A 



Since the evolution of a, (3 and A, decouple from that for 0, we consider the behavior of the reduced 
system 

4a* (3* 2(N c F r D +al) 

-2(N£P + 3a 2 - Al) -4a*A, 

-2A,(3, -2a,A, -2a,(3, 





5a 




5a 




(3.30) dt 


5(3 


= B(t) 


5(3 






. 5A . 











5a 




5(3 







Then, B is the 3 by 3 block in matrix of equation (|3.28|) . In obtaining (|3.30|) . we used that N cr = Qio/2. 
Linearized dynamics about the symmetric equilibrium state: 

For the symmetric equilibrium, as displayed in (|3.22[) we have 

(3.31) (a e _ q ,(3 e _ q 7 A e _ q ,6 e _ q (t)) = (0, 0, JV' , (-O + N)t). 

Hence, 

N[ r D " 
N-NgP 






B = B- = 2 



For the reduced system, we have 
(3.32) B_ = 2 



N FD 

N - N FD 
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Figure 5. Symmetric state for N — N^ r D < corresponds to an equilibrium elliptic 
point (a, (3) — (0,0). Plotted phase portrait corresponds to parameter values N cr = .1 
and N - N[ r D = -.05 (left) and N cr = .1 and N — N[ r D = .05 (right). It is clear for 
N — N^ D < the equilibrium sotion is stable and for N — N^ D > the equilibrium sotion 
is unstable. 



whose eigenvalues and implied linear stability character is as follows: 

(3.33) N < N^ D : A = 0, X±(N) = ±i2 [ (N^ D - N) N^ D ] * , N < N[ r D (stable elliptic point), 

(3.34) N > N[ r D : A = 0, \±{N) = ±2 [ (N - N[ r D ) , N > N™, (unstable saddle). 

see Fig. [5] Thus, the symmetric state transitions from stable to unstable as N increases beyond Nf r D . 
Furthermore, £?_ can be diagonalized 





r 2n%p 


2N™ 


" 




" A+ 





" 




f- 2N™ 


2N™ 







A+ 


X- 










A_ 







A+ 


A_ 













1 





















1 



and the linear evolution is given by 







r 2 N ™ 


2N FD 





(3.35) 




A+ 


A_ 















1 







2N™ 

A+ 




2N FD 

cr 

A_ 




In other words, we have the bound 



8a 



cos(|A + |i)fo 
cos(|A+|t)<5/3- 



\N-N£P\ 

\N-N^ r D \ 
IV™ 

SA 



The full linearized dynamics are governed by the matrix 
(3.36) = 



* * 

* * 
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where B has the same eigenvalues as B, with A = now a generalized eigenvalue of multiplicity two, 
implying linear growth of S9(t). Explicitly, we have 

" N FD (N — N^ r D ) " 

N[ r D (N-N FD ) 








Hence, it is clear 

- B - t - M + imtkh + cos(2y/ N£ D \N - NgP\t)M 2 + sm(2^N^ r D \N - NgP\t)M 3 , 



e 



where 



Mo 



Mi 



M, 



Mo = 







10 

1 







10 

10 

10 







y/\N-N£P\ 






y/\N-N£P\ 


















Note here the leading order behavior is then a system which oscillates at the correct period and grows 
linearly only in the phase term. On a component basis, wc may say 



Alternatively, we can use the spectrum of B_ to define the matrix 

1 10 





1 
1 



Then, we have 



-2i ((N c F r D - N)N, 



FD\ 2 





2i({N™-N)NZ D ) i 











27Vi 
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Hence, we see 



-,-B-t 



P_ 



e 2i((ivf--iV)ivf r °)i t 



1 
1 



p- 



2N?t 
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Linearized dynamics about asymmetric equilibrium states, N > N^ r D : 
For either asymmetric equilibrium, displayed in (|3.22|) we have 



(3.37) {o%,0?,A'«,0°*(t)) 
Hence. 



N-N. 



FD \ 2 



.0, 



N + N, 



FD\2 



(-2f2 + N£ r - N)t 



(3.38) B+ 





-2(N-N™) 




2( 2 




FD\3t r (N +N* r ) j 




^ ( 7V-A^)i(H±£^ 


2(H±^)i 







In this case, we substitute (|3.25|) into the expression for B in (|3.30|) and obtain: 
N + NlP 



(3.39) 



B = 



N + N[ r D 

-2(N - N[ r D ) 2 (N 2 - {N^ r D ) 2 Y 

-(N 2 -(N[ r D ) 2 y 



N > N™ = 



n 



in 



The eigenvalues of B are: 

(3.40) A = 0, \±(N) = ± 2i {N 2 -{N^ r D ) 2 y\ N > . 

Therefore, the bifurcating asymmetric states are stable elliptic points. Also, B + is diagonalizable, resulting 
in 

i(l + co S (|A + |i))fo + (£±$g) h sin(|A + |t)£/3 + ± (f^) * (1 - cos(|A + |t))<L4 



cos(|A+|i)5/3- 



\N-N^ D \ \ 2 
N+N£ 



Sit 1 sin(|A+|i)<5a + sin(|A+|t)<L4 



| ( 1 n+n}P ) ' (! - cos(|A + |i))<5a + ± sin(|A+|t)<5/3 + { + cos(|A + |t))M 
In a manner analogous to the case of symmetric bound states, we have 



' " 









= 0, B+ 





1 





N+N™ 

n-n£p 

1 
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e B+t ff 



r/+sin(2 (N 2 -N 2 cl ) 2 t) (N 2 



jFD N 



(N 2 - N 2 .) U))[ n _ n C fd ) %ei + (N + N™)^! + t) m e 4 . 



Once again, we see the dominant behavior be oscillations in a, (3, A with a linear growth accompanied by 
a factor of order N? in 9. 

As before, we can use the spectrum of B+ to define the matrix 



P+ 



-2i 



1 

N+N* T U 
N-N 



2 1 



1 

N-N* D 

FD 







N+N™ 









and 



PI 



1 / N+N'" 
' i [ N-N*i D 



N+N'' r u 
N-Nf r D 






-J( 

-J I 

i( 



A f +A r cr 
N — N^p 

N+N[ r D 
N-N™ 

N+N™ 
N-N£P 





Then, we have 



P+ 1 B+P+ 



-2i ((N 2 -(N c F r D ) 2 )y 










2i ((N 



T FD\2 



(N c 
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Hence, we see 



„B + t 



P+ 



2i ( ( N™ + N) (N~N[ r D ) ) 3 1 









e -2i{(N^ T D -N)(N-N^ T 









+4x/2(7V - iV c ^)5( e 2 ^ 2 -^)** - i)P + 



PI 



PI 



+4v / 2(^V-^)'(e- 



+8y/2(N-N™)h 



-2i[N 2 -Nl r ) 2 ( 



DP+ 



PI 



2 y/N-N^P 



Remark 3.2. The analysis of e B± * provides partial understanding of the inherent instability in the phase 
of the finite dimensional periodic orbits. Specifically, an orbit (Ai,ai,/3i) will oscillate with period Xi 
while a different orbit (A2, o^Zfe) will oscillate with period T%. As a result, if two orbits begin quite close 
in phase, they will naturally oscillate out of phase with one another, which from the analysis above gives 
the linear shift in the 69 component. However, the oscillations will be purely described by the (A, a, (3) 
system, meaning we must decouple the phase from the equation to get stability. 

We have can the following 
Proposition 3.1. The system (|3 . 1 1[) has periodic orbit solutions. 

Proof. The proof follows from standard level set techniques. Namely, we have the convex geometric curve 
H, which we slice at a particular value of N. Since the a, (3 coordinates are constrained to live on that level 
set, the only possible orbits are those that orbit with a period determined by N. Note, the equilibrium 
solution is the minimum of H with respect to N. 

First of all, we have 



Ha 


= 4a [2a 2 


-(N- 


Hp 


= 2/?plio 


+ 2a 2 ], 




= 4[(6a 2 - 


-(N- 


H (i .p 


= 2[fi 10 + 


2a 2 ], 


Ha,/3 


= 8a/3. 





2 ' 



~2~ 



a 2 



Hence, for ./V < N^ , we have an equilibrium point only at (a, (3) 
equilibrium points at (a, (3) 



(0,0), (± 



N — - 



-,o). 



minimum for H when ./V < N^ r D , and (0, 0) is a local maxima while (± 
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(0,0). For N > N[ r D , we have 
Using the Hessian, we see that (0,0) is a local 
, 0) are local minima when 



N- 
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N < N F r D . Hence, for a fixed N, we select a bounded, closed curve for (a,/3). By the Poincare-Bendixson 
Theorem, the possible curves are asymptotic to a limit cycle. However, the orbit is a closed curve, so it is 
necessarily periodic. □ 

Remark 3.3. We note that the dichotomy in the types of equilibria found for N > Nf r D and N < N^ r D 
exactly describe the phenomenon of symmetry breaking for the profile of the stable bound state for 
at small amplitudes as seen in |21j . 

3.3. Periodic Orbits Near the Equilibrium Point. In this section we estimate the periodic orbits and 
describe them near the equilibrium point. 

Linearizing about the equilibrium solutions of (|3.6|) for N — N FD > 0, one may define 



\N -N FD \ 

(3.41) a = ^ 2 cr l - +h u 

(3.42) (3 = h 2 , 

where hi, h 2 small perturbations. Plugging this ansatz into (|3.6p gives 

(3.43) hi = -4(iV - N FD ) (n fd + hl + 0(\hf 



2 

Hence, the period of oscillation near the equilibrium point is of the form 

(3.44) T = ^ + 0{\N N FD \ 2 ) = T + (N - N FD ) + 0(\N - N FD \ 2 ). 

\1\ +*cr I 

Similarly, for N - N FD < we have 

(3.45) a = hi, 

(3.46) = h2, 
which when plugged into ()6.1() gives 

(3.47) hi = 4|7V - N FD \N FD hi + 0(\h\ 2 ). 
Hence, the period of oscillation near the equilibrium point is of the form 

(3.48) T= r + 0(\N - N FD \ 2 ) = T-{N - N FD ) + (D(\N - N FD \ 2 ). 

V I Q N _ N FD\ N FD^ Vl c r I ) V cr ) VI 

It is on a time scale of multiple oscillations we hope to control the difference between the observed finite 
dimensional periodic solutions and the full solution to the PDE. 

This leads us to the following: 

Proposition 3.2. Fix (N - N FD ) > or (N - N FD ) < such that \N - N FD \ <C 1. Let us assume 
N F r D 3> \N — N FD \. Take (Po 9 , Pi 9 ) be the corresponding equilibrium solution for (|3.1|) . For any periodic 
solution po(t), pi(t) of (|3.ip such that 

(3-49) IO*(*),Pi(*)) - (Po q ,Pi q )\ « \N~N FD \, 

we have 

(3.50) (p (t + T),pi(t + T)) = (p Q (t),pi(t)), 
where 

(3.51) \T -T±(N- N FD )\ « \N - N FD \ 2 . 
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Remark 3.4. For small perturbations of the equilibrium point (for either (N — N^ r D ) < or (N — N^ r D ) > 
0, it is for precisely the period T±(N — N^ r D ) on which we must control the coupling to the continuous 
spectrum for the full solution to in order to prove these finite dimensional structures are observable 
over many oscillations, see Section [5] In order to generalize our result to any periodic solution predicted 
by the finite dimensional dynamics, we must understand fully the period of each full oscillation. This will 
be discussed further in Section [7] 

Finally, we can state the following 

Proposition 3.3. Fix e > 0. There exists a S > such that if a given periodic orbit solution of (|3.26p . 
(A(t), a(t), (3(t), 0{t)), with period T with 



(3.52) | (,4(0), a(0), 0(0), 0) - {A eq , a eq , f3 eq , 9 eq )\ < 6, 
then we have 

(3.53) \\M(t)\\L-^L- < (l + ^He^lU-^L-, 

(3.54) \\M(t)\\ L oo^ L oo < (l + eJUe^lli-^i-. 



Proof. The proof follows from continuity of the Floquct multipliers. We discuss the case for B eq and M(t) 
here. The analysis for the case for B eq and M will follow similarly. 

It is clear that (A(t),a(t), $(t), 0(t)) is a solution to (|3.28p . giving at least one Floquet multiplier Ai = 1. 
Similarly, diffentiation with respect to the period gives a similar result, meaning we in fact have Ax = A2 = 1. 

From the analysis of the phase diagram for periodic orbits in a, 0, we have that a(t) = a(—t) and 
f3(t) = -/?(-*). Hence 

(3.55) / a(s)P{s)ds = 0. 

Jo 

So, given tv(B(t)), we know from standard Floquet theory 

4 

(3.56) [] X J = L 

3=1 

Hence, either A3 = A4 with | A3 1 = | A4 1 = 1 or A3, A4 G R with A3 = AJ 1 . By continuity of the Floquet 
multipliers, near cS eg , we have the only the degeneracy at A = 1 resulting in a similar growth behavior to 
that of e B " qt as seen in Section [XU □ 

4. An Ansatz for the Coupled System 

Starting with the finite dimensional system given in (|3.7|) . we consider a periodic orbit, (a(t), /3(f)), near 
equilibrium point and construct a periodic orbit of the extended system (|3.6[) : 

Below we shall specify how near the equilibrium <t» need be. 

We now write the system (|2.22[) by centering around the orbit cr* of the finite dimensional truncation: 

<7(i)=<7*(t)+77, 

(4.1) = (A(t) + ?M (t), a{t) + Va {t)J(t) + T) (tfj . 
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This corresponds to a solution of the form: 

u(x, t) = e ie ^ + rj A (t)) Vo + [(5(t) + »/„(t)) + i(fi(t) + r]p{t 

with initial conditions 

u (x) = e i6 ^ ( i(0) Vo(x) + [5(0) + 0(0)] V'i(z) 

Centered about ct* , the system (|2.22[) becomes: 

'7 = D a F FD (<r*(t)) rf+ [F FD (a* + rf) - F FD (a*) - D 3 F FD (a*)rf\ + G FD {a*,rf; R, R), 
6 = -fi + A 2 + 3a 2 + P 2 + Ge{R, R; A, a, (3), 
iRt = (H- n )R + {A 2 + Za 2 + /3 2 )i? + F b (a, + r,) + F R (a* + ij; R, R) 



+ IA—A + 3(a z -a z ) + [3 2 - (3 Z J R . 

Wc have the estimates: 

Ffd{5* + ff) - Fpjj (<7*) - D ff F FD (^)rf | = ( i^ 2 + |r/| 3 
See Appendix [X] for the explicit expressions in the system. 

Wc denote by R the leading order part of R, driven by the periodic solution a*(t): 

(4.2) iR t = (H- n Q )R + (A 2 + 3a 2 + p 2 )R + P c F b (a„). 
The correction to R is given by w, which satisfies: 

(4.3) R = R + w. 

Introduce, M(t), a fundamental solution matrix for the system of ODEs with time-periodic coefficients: 

dtv = D a F(a«{t)) r). 



We shall study the following system of integral equations for T](t), 9(t), w(x, t): 



ff(t)= / M(t)M- L {s) (F FD {a* + rj) - F FD (a*) - D a F FD {a*)ff) + G FD (a*,m R, R) 



ds, 



(4.4) 6 = 6»o+ / [-n + A 2 +3a 2 +f3 2 +G e (R,R;a*,ff)]ds, 
Jo 

w(t) = f\imt-s)-ia (t- s )+if^ +3Sl ^)( s ') ds ' Pc [ {Fb &+ff)-F b (^)) + F R (a*,ff;R,R)] 



We view a solution, (if, i?), of this system of integral equations as fixed point of a mapping, M: 
(4.5) (ff,R)=M(ff,R). 

In the next section we formulate and solve this fixed point problem in a function space, which yields 
the existence of solutions to NLS/GP which "shadow" the periodic orbit, <r*(£), on the time scale of many 
periods. 
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5. Main Results 

Recall that the period of the orbits described in Proposition 13.11 satisfies 

T penod ~ (\N[ r D -N\N c F r D y k i. 

In this section we shall construct solutions to the full PDEs, which "shadow" the finite dimensional orbits 
for many periods. 

Let r > denote a number to be chosen sufficiently small. And define the region of parameter space, 
about the symmetry breaking point, in which we conduct our study by 

\N FD -N\ = t, 

(5.1) N c F r D = t\ 

For now < 7 < 1, but we shall place further constraints will be placed on 7. 

Remark 5.1. Recall that r = N FD ~ Af cr (L), where L is the well-spacing parameter for the double-well. 
Since, for L sufficiently large, N cr {L) ~ &i(L) — Qq{L) < e~ KL , the eigenvalues splitting , and r can be 
made small by choosing L sufficiently large. 

In terms of r, we now describe the periodic solutions discussed in Propositions 13. ll and [3~2"l 



Proposition 5.1. There exists tq > 0, such that for < r < tq, the system ()3.7j) has periodic solutions, 
which are small perturbations of the equilibria: 

(5.2) \Mt)\ 2 -r\ \a{t)\ 2 + \P{t)\ 2 ~ r, N > N c F r D , 

(5.3) \Mt)?~T\ \c\{t)\ 2 + \~P{t)\ 2 ~r 1+s , N<N FD , 
where 6 > is to be chosen below. 

The period of oscillations of the periodic solutions of Provosition YS '. 1\ is 

(5.4) T penod (T) - t- 1 ^. 

We shall seek to contruct solutions on a time interval of the form 

(5.5) 0<t<T*(r), T^t) ~ r~ e x T pe „ od (r) = r^" 6 . 

5.1. Notation. It should be noted, notationally when we refer to 

(5.6) A<B, 
we mean 

(5.7) A < CB 
for some C of order 1. Also, for 

(5.8) A<B, 
we mean 

(5.9) A < cB 
for some c > much less than 1. 
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{a eq ,p ev ) 



Figure 6. Let N > N^ r D . By Theorem Q] for data (a ,/3o) on the periodic orbit 
( a(t),(3(t) ) and such that |5o — ot eq \ 2 + \fio\ 2 ~ t 1+s , the solution u(x,t) evolves within 
the dashed annulus of width t 1+Si , 5i > S. 



5.2. Statement of Theorem. 

Theorem 1. Assume < r < tq and | < 7 < 1. 

Denote by 



a periodic solution of (|3.6[) . /or which: 

(1) whose period T per i 0( [{T) satisfies 

T pmod (r) < (\N c F r D - JV^™)"* = r- 1 * 2 

(2) whose fundamental matrix, M(t) of the linearized dynamics about o~*(t) satisfies the norm bound: 

( N + N FD V 7-1 

= C 



0<s,t< T perlod (r) ||M(t)M- A («)|| < C ^ 

(3) and which is a small orbit about an equilibrium ( A eq , a eq , (3 eq ) as in Provosition \5 . 1\ See figures 
H and 

Tafee initial data for NLS/ GP of the form 

(5.10) uo(x) = e^ ' ( i(0)Vto(x) + [5(0) + 4/9(0)]^ (*) ) , 

where 9(0) G R is chosen arbitrarily. 

Then there exists a solution u(x,t) of 0/ f/ie /orm 

u(x, t) = e M <*> + VA {t))Mx) + [(5(*) + »/«(*)) + W) + (a) + t) + w(x, i)) , 

where R as in (J4.2I) . 

Furthermore, T)(t) = ( ■q A {t),n a {t),i lp (t) ), 0(i) G ^([O, T*(r)]), to G L^H^ DLfLf and w satisfies the 
bounds 

(5-11) IM|Lf[0,T„(r)] + lk|U~([0,T»(r)];ffi) + || w\\ L 4([ ,T. (t)];£~) ^ , 

/or all t € I = [0,T*(r)] = [0, T per i OC i(T) r _e ] where e > is specified in the proof. 
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(a eq ,f3eq) 



Figure 7. Let N < N^ r D . By Theorem Q] for data (a ,A)) on tnc periodic orbit 
( a(t),/3(t) ) and such that |<So| 2 + \Po\ 2 ~ t 1+s , the solution u(x,t) evolves within the 
dashed annulus of width t 1+Si , Si > S. 

5.3. Proof of theorem. In the sequel, we will often use the contracted notation 

L\Wl' q = L p ([0,T*};W k - q (R)), 

where T* is given by (J33J), (f5TT|) . 
We now show that the map 

(rf,w) = M(ff,w) 

defined in section [4] is a contraction. 
Define the space 

X(I)=X([0,T*(T))) = {(n,w):r,€L™([0,T*(T)]), w e Lf ([0, T,(r)]; i£) n L*([0, T*(t)]; L~) } 
equipped with the natural norm 

ll(»7»IU(/) = II»7||l?°(7) + I|w||l*(7;l-) + l|w|U r (/; ff i), 

where /=[0,T»(r)]. 
We define B T (I) C X(J~) such that (77, i?) G B T (I) if and only if 

(5.12) \\{ri,w)\\x { i)<T-i +S \ 
where 5\ > will be chosen later. 
We must prove the following: 

Proposition 5.2. T/ie mapping Ai : X(I) — ► X(I), defined in (|4. 5[) . /ias i/ie properties 

(1) X : B T (J) -» B r (I). 

(2) There exists k < 1 skc/i t/iai given (ffj,uij) G B T {I) for j = 1,2, we have 

d(M(ffi,w 1 ),M(ff 2 ,w 2 )) < k d{{ff 1 ,wi), (if 2 ,w 2 )). 
Thus, there exists a unique solution (jf, w) in B T (I) . 

Proof. We begin by proving necessary bounds on R. 

24 



J.L. MARZUOLA AND M.I. WEINSTEIN NLS / GP WITH DOUBLE WELL POTENTIALS 

Proposition 5.3. Given a solution <?» = (A(t), a(t), f3(t)) to (|3 . 7[) . as in Proposition \5.1l there exists 
To > 0, Sq > S, such that for all r < tq we have the following: 
for all t £ / , and defining 

(5.13) R= e -iH(t- s )-ia (t- s )+if^A 2 +3& 2 +^ 2 )( s ')ds' PbF ^ )ds 

Jo 

for Ff, defined in (|2.1[) , we have 



\\R\\l?>(I;L«>) 



< T 1+s " 



Proof. Note that it is only the magnitude of the components of <t*, via F(<j*(t)), which factor in to the 
bounds of R. 



From Proposition 15.11 and the finite dimensional conservation laws 

(5.14) A(t)=rl + eo(t), 

where |eo| ~ t~%~ . Based on the expansion in the ansatz, we have 

^0*) = Pc[A\>$ + (i 2 (<5 - i[3) + 2A 2 (& + *0))VoVi + ((« + + 2i(a 2 + ~P 2 ))Ml 

+ (a 2 + f3 2 )(a + if3)ipf]. 

The term of largest order in this expansion is 

(5.15) A 3 iPl 

The bounds on the remaining terms will follow similarly, so we look at 



-i(H-no)(t-s)-iSZ(A 2 +3a 2 +p 2 )(s')ds'p c ( T % + eQ )3^3 ds 
_ T ^ f e -i(H-n o )(t-a)-if^A 2 +3& 2 +0 2 ){s')ds' p c ^3 ds _ 

Jo 

In particular, we will show there exists So > 5, such that 

(5.16) ||r* [\-« H -WM- i M A *^+PWWp e $d8\\ L - <t 1+s °. 

Jo 

From p.6p . we know 

(5.17) §= -tt a + A 2 + M 2 (t) + [3 2 . 

Hence, the leading order constant terms from the derivative are N^ r D — Qq. We write 

e -i{H-n ){t- S )-tfl(A 2 +3a 2 +f3 2 )(s')ds'p c ^3 ds 



T 







e -iHt-i8(t) I e i9{s)+in s-iN^ D s e iHs-in s+iN^ D s p^3 ds 



ds i 



o 
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where the resolvent (H — Slo + N^P) 1 is a well-defined operator in H 1 on P c ipo, see Appendix IBl Using 
the estimates from Appendix [Bl and integration by parts, we have 

-iHt-iG{t) f e i0(s)+in o s-iN[ r D s e iHs-in s+iN? r D s p^fa 



II 



2± 



+t^O(t) f \\e iHs (H - fi + N^y 1 P c ^\\ L ~ds 
Jo 

< t*\\(H Q„ + N c F r D y'P c Mm + r^\\e- im e-^{H - fi + N^P^Wm 

+t*0(t) / \\e lHs (H -Qo + N^y'P^IW^ds 
Jo 

3 1 37 3-Y , . 

< T~ +T~ + T~0(t) t. 

By selecting 7 > |, we ensure that all terms resulting from integration by parts are bounded by t 1+Si for 
alH G I = [0,t-^- £ ]. □ 

Using dispersive estimates, we have additional bounds for R given by the following 
Lemma 5.4. Given R defined as in (|4.2[) . we have 

\\R\\m<\A\ 3 \I\i 

and 

\\R\\L H i;Li) < 

for any Strichartz pair (p, q) . 

Proof. To begin, we first note that by the L 2 boundedncss of wave operators discussed in Appendix [Bj we 
need only consider 

\\(H)R\\li- 

As in Lemma 15.31 we look at the the worst term, namely 



e ~iHt e -i$(t) 



(H-n Q + N^y'p^ 



(H S )R\\ L ~ L 2 < 



e *B(t-')A 3 (a){H")P c il4da\ 



L^L 2 



< T*\I\*, 

where we have used the Strichartz estimate (|B.10[) with dual Strichartz norm L%(I; L 1 ). 
The Strichartz estimate follows similarly. 



□ 



Remark 5.2. Though such estimates do not arise in this work, let us also point out the following simple 
estimate 

r-i 



(5.18) 
(5.19) 



\(R,Xi}\ 



< 



3j_ 

T 2 







< T 2 



{t — s) 2 ds < t 2 t 2 , 
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which is proved in |32| and may be applicable when trying to prove long or infinite time results on similar 
problems to the one studied here. 

Note, by standard Sobolev embeddings and the contraction assumption, we have 
(5-20) p?|U~ < T i+*i. 

Since we have proper bounds on R, we must now bound 

IMIl- = \\ f M(t)M-\s)[F FD (<7,+rf)-F FD (a,)-D a ,F FD (a,)rf 
Jo 

+G FD (a* +rf;R,R)]ds\\ L <~ 

and 



t 

iH(t-s) jn (t-s) i J*(A 2 +3& 2 +0 2 )(s')ds' 



\\w\\L<x>HinL*L°° = !l / e v 'e >e 
Jo 

xP c [(F b (3> +rf) - F 6 (ct*)) + F R (a*,rj; R, R)] ds\\ L ^ H i nL i L ^. 
From Appendix [Bj we have for any Strichartz pair (p, q) that 

(5.21) || f e iH ^PJ\\ LPW i, 9 < \\f{x i t)\\iSwi<*> 

Jo 

where (p, q) is a dual Strichartz pair. In one dimension, it is useful to note we may take p = | + \i\ and 



q = 1 + fj,2 for [ij > small, j = 1,2. In other words, we can be as close to the endpoint estimate of p = | 
and q = 1 as we need to be. 

By reducing the system and including the lowest order terms from the expansions above, we can reduce 
to controlling a model problem of the form 

* = [ t ( N+N ZXww 

1 J \\N-N£P\) [ 1 1 

+ A 2 (R, X ) + A 2 (w, X ) + 2MX, R) + V 2 (R, X) 

+ At}{w,x) +V 2 (w,x)]ds 

= Ternil + Terrn^ + ■ ■ ■ + Terrrig], 



w 



where and 

t 

e iff(t-s) e m (ts) e « f*(A 2 +3& 2 +0 2 )(s')ds' 

x P c [(3i 2 ??Vo + ZARt] + Aip 2 Rr]) 

+ {A 2 ^lw + w 3 + 2Ariw + A(x, w)R + A(x, R)w) 

+ (A 2 ij 2 R + A( x ,R)R)}d s 

= 1 erm 1 + 1 erm 2 +■■■ + ! erm 10 . 

For our model problem, we now take (77,1/7) G X, \\(rj,w)\\x < t^ +Si . In addition, take (i]j,Wj) G X, 
|(^-,^)IU < ^ +Sl forj = 1,2. 
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It follows 

/ N + N FD \ * /"* - 

< ^ 2 + 2 2 e Nl£« 

provided we choose 

(5.22) 5l> L-l + e . 

In addition, it is clear from the analysis there that 

|| TermY ~ Term 1 ] 2 1| l=o < hl5l_c ||^i — 772 1 1 i = 

Next, for Term\ we get 

/ at 1 /uFD \ 5 /■* 

» Te ™"»- - (jlv^l#|) II / 'W*"" 

7 1 1 + 7 .. Q 



< T 



H 3,5 i- e 



provided 

(5.23) 2<5i > e. 

Similarly, we have 



\\Termf -Term? < ^r 25l ~ e ||r7i - ?? 2 |U= 



For Terml, we have 



llrermglU- = ( |jV „^| j II ^ ^ 2 <^X>^IIl= 

provided 

(5.24) 7 _ e> I +( 5 1 . 

Since Terml is independent of (77, w), it follows easily that 

HTermg 1 - Terrr^ 2 j| L oo = 0, 
meaning this term does not factor into the contraction. 
The bound on Term\ is of the form 

/ N + N FD \ ^ /"* 

ll Te ™"»- = {w=W0 I J, 



< 7-2 8 4 
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provided 

(5.25) 1> \ + \ t - 



It follows immediately that 

\\Termf -Termf || L « < r^+T-f-^-T || w || i4L o 



< T S V _ W 2 j| L 4 L o 



For Terrnl, we have 

/ N + N FD \ * /"* - 

< T*-*+i-*-i- e |WU-||iZ||i- 

< ri- e h|| L oo < T 5+ 51 , ( 7 >2e), 
which follows from previous constraints on 7. Once again, we have as well 

\\Termf -Termf\\ L oo < ri^\\ m - m \\ L ~. 

For Terrnl, we have using Proposition [53] 

/ N + N FD \ ' /"* o - 
HTermglU- - ( |jvljv|^| j 11 /„ 'ft*'*)* 

< r i-Hi-l-i-||, 7 ||| 00 

< r 3 +5l - e ||?7||l,=o < T $ +Sl , 

which follows if 5i > e. Furthermore, 

\\Termf - Termf \\ L ~ < ri +S ^\\ m - 77a || L -. 



For Terrn\, we have 



[ \ n _ n fd\ ) II / Mw,x)ds\\i 

< 2 + 2 8 8 4 |j7 ? |j L30 |j W || i4L = 

< r ^-|-^+5 1 || r? || Loo < r i+*i 

provided 

3 6 8 r 

(5.26) 7 > - + -£- -5,. 



Furthermore, 



5 5 5 



1 1 Term? 1 -Term? \\ L °° < r^""^ || ( m - m , Wl - w 2 )\\ x . 
For Terrnl, we have 

jFD \ 5 ft 



< T2 2 S 8 4 ||7?||2,ooH|z,4I,°o 
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by previous constraints on 7. Again, it follows that 

WTerml 1 -Term v s 2 \\ L ~ < r*+*-* +2<51 |K - w 2 \\ l a l ~>. 



We now study the map on the dispersive part, w. For Term™ , using the Strichartz estimates (|B.10[) , we 
have 

HTerm? Ui-^n^cc = || f jW->)jW-*jr. a>+**+P)Wd. ' p c3 A 2 r,^ ds\\^ mnL ^ 

< ^-i-^U\ L ^<T^ 5 \ ( 7 >jj), 
by previous constraints on 7. Hence, 

\\TermT -TermT\\ L ~mnL*L~ < t^"*"*-* H771 - tja|U~- 
For Term£ , we have using bounds on i? in if 1 (Lemma [57 



HTermyiliootfin^i- = || / e lff ( t -^e j0 °( t -^e I ^( i2+352+ ^^ s ') ds 'P c 2iii7 7( is|| 



< T*\\R\\ L i H i\\r,\\ L 

< 

which follows using 7 > | + -^e. Hence, it follows that 



7 77 7e .. 97 7 7e ., 1 , c 

< r 2 T-s-T--|| ?7 || icc < r--B--||r y || L =c < ,-*+*, 



\\Terrri2 1 - Term% 2 \\ L ^ H i nL i L ^ < t * I » \\r)i - ?/ 2 ||l~- 
For Term™ , using Lemma 15.41 once again we have 

\\Term w \\ L - HWL ~ = \\ f* e^.^iOo^i^JP^a'+^d. ' p c A^Rr,ds\\ L ~ HWL ~ 

Jo 

< t1\\R^\\ 

< T^-i-^-^ML- < r T?-S-*||,,|| ioo < T ^ +s \ 
which follows easily from previous constraints on 7. Hence, it follows directly that 

\\Term Wl -Term W2 \\ L ^ H i nL A L ^ < r^~J~^ \\r)i - 7/2||i«. 
For Term w , we have 

\\Termy\\ L ~ H , nL . L ~ = \\ f e mt-s) e iU (t-s) e if:^ + 3a^)( s - )ds ' p c A^ oW ds\\ L ~ HWL ~ 

Jo 



which follows from previous constraints on 7. Furthermore, 



\Term™ 1 - Term™ 2 \\ L ^ H i nL i L ^ < t« + 2 £ \\ Wl - w 2 \\l°°h 1 ■ 
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For Term§ , we have 



ni 4 L° 



— i — i 

< T 2 2 



which follows easily by previous constraints on 7. Moreover, 



For TerrriQ , we have 

||Terr<|U»ifW£» = || /" e ^-') e ^i*-') e *St ( A2+3&2 +^^ ds ' 2P c Ar 1 wds\\ L ^ m 

Jq 



ni 4 i° 



which follows easily from previous constraints on Si. Hence, it follows directly that 
WTerm^ 1 -Term,g 2 \\ L oa H i nL 4 L oa < T Sl " f -\\wi - w 2 \\ L ^H 1 ■ 

For Term™ , we have from Lemma 15.41 



nL 4 L° 



_j _! j_ 



< 2 2 e \\R\\ LrH i\\ w \\ L °, H i 

< T^-i-^Hwll^^i <T^', (7>g). 
which follows from previous constraints on 7. It follows directly that 

9-C 7 7 

HTerm™ 1 - Term™ 2 \\ L ^ H i nL i L ^ < T^~s~i e ||u>i _ w 2 || L ^ ff i . 

The bounds for Terrrig follows in a similar manner. 
For Terrrig and Term™ , we have 

Jo 

and 

llTermrolU-^n^oo = || /'^M^lf-.)^ (2 2 + 3 a = + ^)( s ')^ 'pj^ fl)Afa|| iooHlni4ioo 

Jo 

< T-THfllliifl-illfllli- < r 27+1 -*T-*-3 e < T5+ 5 '. 

Both of these terms are independent of (rj,w), meaning the contraction mapping follows easily. 

Hence, choosing 7, Si, and e such that the constraints (|5.22p . (f5T2"5)) . (|5.24[) . (|5.25p . (|5.26[) arc satisfied, the 
contraction argument follows and the result holds. 

Once we have solved for (a,R), it is clear that the resulting function 9 £ C(I) by construction. □ 
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6. Numerical Simulations 

6.1. Polar Coordinates. As it will simplify the process of building initial conditions for numerically 
solving (jl.lj) that display the behaviors we study above, let us discuss here an alternative set of coordinates 
for (|3.1|) . Namely, we set po = roe 1 ® and p\ = r\e %dl . This leads to the following system of ODE's: 

r r = r?r sin(2A0), 

(6.1) I r\ = -r%n sin(2A0), 

{ (AO) = fii - il Q + (rf - rg)(l + cos(2A0)), 

where A# ~ 0\ — Oq. Given the system above, we can say that the bifurcation of stability occurs at 
(6-2) N™ = 

As we are interested in the behavior quite near the bifurcation point, we define new parameters eo, ei and 
n such that 



(6.3) r = v /iV^+e , 

(6.4) n = ei, 

(6.5) A = aJ; z? + ? 1 , 

where 

(6.6) n = eg + e\ + 2^N^e^. 
Then, wc have 



(6.7) e = 4 [jNF r D + e sin(2A0), 



(6.8) ei = - ei J A™ + e sin(2A0) 



2\ 




(6.9) (A0) = n x - fio + c? - V^™ + eo (1 + cos(2A#)) 



(6.10) 



It should be noted, using the conservation laws we can write the system for ei and AO independently 

e\ = -ei (N™ + n-e\) sin(2Afl), 
(AO) = Oi - Q - ((a£ d + n - 2c?)) (1 + cos(2A0)) 

and hence analyze phase plane diagrams, see Figures [5] and [5J In particular, the behavior of A8 in these 
regions leads to interesting oscillatory behavior as shown in the phase diagrams featured in Figures [5] and 
IH1 For n > 0, a simple calculation shows that the equilibrium solutions occur for e\ = y^f, AO = kir for 
k e Z. 

For n > 0, we have trapped orbits near the values AO = kir for k 6 Z. This is a manifestation of the 
orbital stability of these mixed states. However, oscillations between wells can be generated by a large 
enough phase shift to leave the region where such trapped orbits occur. These oscillations are then large, 
in particular e\ must reach some e max before decreasing. 

For n < 0, we see the oscillations still attain a maximum at AO = kn, however their amplitude approaches 
with the initial ei(0). This is a manifestation of the stability of the symmetric state in this regime. One 
may ask if such finite dimensional Hamiltonian dynamics would appear in the infinite dimension dynamics 
of the PDE. For this, see Figures [TUl [TT1 and [T2l for numerical evidence of their existence for long times. 
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Figure 8. A numerical plot of a phase plane with A9 plotted versus e\ for n => 0. Here, 
we have N^ r D = .2, n = N - NlP = .05. 



Figure 9. A numerical plot of a phase plane with AO plotted versus t\ for n < 0. Here, 
we have N£P = .2, n = N - Nf r D = -.05. 

To begin, we locate the symmetry breaking point for a particular system. To do this, we use spectral 
rcnormalization with an asymmetric initial function to find the symmetry breaking point, sec Fig. [3] 

In the remainder of this section, we run simulations with wells of the form 
(6.11) V L = q(S(x-^)+S(x + ^)). 

Similar results hold and the same numerical analysis tools arc applicable for potentials with more regularity. 
Knowing the bifurcation point as discussed in [17) . we can numerically integrate using a finite clement 
method similar to that in [T3], where scattering of soliton solutions across single delta function potentials 
was analyzed (see [2] for analysis of finite element methods for nonlinear Schrodinger equations without 
potential). It is quite simple to adapt the method presented there to allow for a double- well potential (for 
delta functions or smoother potentials). The initial data is generated by finding the lowest entries of the 
spectrum of the discretized representation of H — — A + V in the Galerkin approximation, which is an 
operation embedded in many numerical software programs. For simplicity, we use the eig function from 
Matlab. Then, we may numerically solve the PDE system ([1.1)1 with initial data corresponding to that 
necessary for the three types of oscillation described in Section [3] Note, one could also use the solitons 
from the spectral renormalization code (see [6]) as initial data quite easily, however these represent true 
nonlinear structures and we wish to observe structures derived from the finite dimensional dynamics, which 
we only expect to persist on finite time scales due to the nonlinear structure. The orbital stability of the 
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Figure 10. At top, a numerical plot of the phase plane diagram for n = N — N^ D > 
with specific points chosen along a closed orbit which shows full oscillation of mass from 
one well to another. Below numerical plots of the absolute value of the solution to Equation 
p.ip at various times with initial data such that n = N — N^ r D > and A8(0) = 1. The 
plots correspond to points a, b, c, d, e, and f respectively from the specified orbit. 



nonlinear objects is an interesting question in its own right and was explored in [21] . The equilibrium point 
of our dynamical system is in fact the finite dimensional part of a soliton solution, so there is no question 
the orbital stability of the soliton and the long time existence of oscillations near an equilibrium point are 
related. The phase plane diagrams for the finite dimensional dynamics are plotted using the MATLAB 
software program pplane7 [3]. 
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Figure 11. At top, a numerical plot of the phase plane diagram for n = N — N^ r D > 
with specific points chosen along a closed orbit which shows localization of the mass on one 
side of the well. Below, numerical plots of the absolute value of the solution to Equation 
(fTTTjl at various times with initial data such that n = N - N^ r D > and A0(O) = 0. The 
plots correspond to points a, b, c, d respectively from the specified orbit. 

Finally, by taking a system such that |f2o — S^i| is comparably large (or well-separation distance, L, 
comparably close but sufficiently large to guarantee the hypothesized discrete spectrum), we can observe 
for large perturbations coupling to the continuous spectrum and hence decay of a fully oscillatory solution 
to a ground state in a finite time. In general, the mass dispersion is rather rapid, hence after a prescribed 
number of time steps determined by the computational domain, we cut off the solution near the origin and 
continue solving with the cut-off initial data. See Figures [131 [TJl [15] and [16] for various computed time 
evolutions of the (radiation) damped oscillations observed in such a system. 

7. Conclusion and discussion 

The fact that we may observe oscillations between potential wells in systems such as (|1.1[) is not a 
new result, however we have been able to give a representation of the phenomenon in terms of a classical 
oscillatory system. In future work on double well potentials, the authors hope to prove long time stability 
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Figure 12. At top, a numerical plot of the phase plane diagram for n = N—N^P < with 
specific points chosen along a closed orbit. Below, numerical plots of the absolute value of 
the solution to Equation (jl.lj) at various times with initial data such that n = N—N^ r D < 
and A6(Q) = 0. The plots correspond to points a, b, c, d, e, and f respectively from the 
specified orbit. 



for oscillations far from the equilibrium point and optimize the time of existence proof by having better 
control of the damping caused by coupling to the continuous spectrum. In addition, these pseudo-bound 
states represent possible solutions resulting from the problem of scattering of solitons across double well 
potential wells, which in the high velocity limit have been recently studied in [1]. 

Finally, the authors would like to point out this question of oscillation and the resulting dynamical 
systems becomes more challenging and interesting as the number of wells is increased, see [20] . In particular, 
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Figure 14. A numerical plot of the location of maximum amplitude of an oscillatory 
solution decaying to a ground state and in particular the transition region blown up to 
capture the fast oscillations. 




Figure 15. A numerical plot of the center of mass of an oscillatory solution decaying to 
a ground state and in particular the transition region blown up to capture the fast oscil- 
lations. 

as one increases the number of wells to oo. the phase shift required to see oscillation from one well to the 
next might give insight into the celebrated Pcireles-Nabbaro barrier for discrete nonlinear Schrodingcr 
systems, see [27], [22], [H]. 

Appendix A. Error Estimates for the Finite Dimensional Ansatz 

We assume we are near the symmetry breaking equilibrium point, meaning we may take a(t), f3(t) <C A(t) 
and A(t) > c > 0. Then, we have 
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Figure 16. A numerical plot of transition across the scparatrix for the full infinite dimen- 
sional problem in both the projection onto the phase diagram and a plot of the amplitude 
with respect to time. 



(iA -OA- n a A)e l9 ijo + (ia-/3- 9{a + i/3) - (a + i/?)^i)e l <Vi + iRt - HR -9R = 

- [A 3 ^ 3 + (a 2 + (3 2 )(a + iftifc + A(a + i/?) 2 V>iV'o+ 
2A(a 2 + 2 )^o + A 2 (a - i(3)^ 1 + 2A 2 (a + iflipfy 1 ] 
- [2A 2 4>1 + 4A(a)Mi + 2(a 2 + f3 2 )^' 2 ] R + [A 2 ijj 2 + (a + i(3) 2 ^ + 2A(a + i/3)ip ipi] R 

- [Ai> + (a- i/3)^} R 2 - [2Aip + 2(a + i/3)^} \R\ 2 - \R\ 2 R. 

Hence, 

(iA — OA — n Q A + A 3 + A(a + i/3) 2 + 2A(a 2 + (3 2 )) = 
- [A(iJj 2 , R 2 ) + {a- i/3)(Mi,R 2 )] ~ [2A{1%, \R\ 2 ) + 2(a + i{3)(Mi, \R\ 2 )] 

- [(%,\R\ 2 R)] . 

As a result, we have 

(ia-/3- 9(a + i[3) - Q, 1 (a + i(3) + (a 2 + P 2 ){a + i/3) + 2A 2 (a + i/3) + A 2 (a - = 
- [A^o^R 2 ) + (a- 1(3)^1, R 2 )] e i2f> - [2A(^ 1: \R\ 2 ) + 2(a + i/3)(^ 2 , \R\ 2 )] 

- [^o, \R\ 2 R)] 

and 

iRt -HR- OR = -P c [A 3 ^l + (a 2 + /3 2 )(a + 0)^1 + A(a + i(3) 2 ^ 
+2A(a 2 + /3 2 )yj 2 ij a + A 2 (a - ifi)^ 1 + 2A 2 (a + i(3)^] 
-P c [2A 2 i> 2 + 4A(a)^i + 2(a 2 + [3 2 )^ 2 } R- P c [A 2 ^ 2 + (a + if3) 2 ^ 2 + 2A(a + if3)Mi] R 
-P c [Ai, Q + (a - i/?)V>i] R 2 - P c Wo + 2(a + i/3)Vi] |i?| 2 - P C \R\ 2 R 

or 

iRt -HR-0R= [F b (A, a, (3, 0) + F R (A, a, (3, 9; R, R)] , 
where we have assumed 

P c F b = F b , P C F R = F R . 
38 



J.L. MARZUOLA AND M.I. WEINSTEIN 



NLS / GP WITH DOUBLE WELL POTENTIALS 



Let us take A > 0. Then, we see 

A = -2a/L4 + Error^, 

a = [Tli - {a 2 + /3 2 ) - A 2 + 9]P + Error^, 

= — [fii — (a 2 + [3 2 + A 2 ) — 2A 2 + 9]a + Error^, 

AO = -Q A + A 3 + A(3a 2 + 1 ) + Error' e 

iR t - HR -6R = F b (A,a,P) + F R (A,a,/3;R,R). 

Specifically, we have (|2TTT]) - (f2TT5|) with 

Errors (R, R, a) = Im (-([2A 2 i<l + 4A(a)VoVi + 2(a 2 + /3 2 )V?V>o] , #) 

- ( [A 2 ^ + (a + »/?) V 2 ^o + 2A(a + i/?)V>3Vi] , i?) 

- [A^lR 2 ) + (a-i{i)(^ u R 2 )} 

- [2A(^ 2 0l \R\ 2 ) + 2(a + i/3)(^oVi, l^?| 2 )] + (1>o, \R\ 2 R)) , 



and 



Error Q (i?,i?,a) = Im (-([2AVoV'i + 4A(a)ip ipl + 2{a 2 + /3 2 )^] , R) 

- ( [A 2 ^ + (a + iP) 2 i>\ + 2A(a + i(3)Ml] , R) 

- [A(M u R 2 ) + (a-if3)(^i,R 2 )] 

- [2A(^i, \R\ 2 ) + 2(a + z/3)(V 2 , \R{ 2 )] + fa, |i?| 2 #)) 

- /3^ 1 Re(-([2A 2 ^ + 4A(a)^ 2 ^i+2(a 2 +/3 2 )V' 2 V'o] ,-R) 

- ( [A 2 tJjq + (a + z/3) 2 ^ 2 + 2A(a + i0)^] , R) 

- [A^lR^ + ia-iPMo^R 2 )] 

- [2A{^ 2 , |i?| 2 ) + 2(a + i/3)(^oVi, l^| 2 >] + ftfo, |i?| 2 i?» , 

Error^JJ, i?, a) = - Rc (-([2A 2 ^Vi + 4A(a)VoV> 2 + 2(a 2 + /3 2 )V>?] , i?) 

- ( [AV 2 Vi + (a + *j9) Vi + 2A(a + l /?)V'oV' 2 ] , R) 

- [A(iPoipi,R 2 ) + (a-if3)(iPlR 2 )] 

- [2A(i> Tpi, \R\ 2 ) + 2(a + i/?)(^ 2 , |i?| 2 )] + (Vi, |i?| 2 i?)) 

- aA- 1 Rc(-([2A 2 ^ +4A{a)^ 1 +2{a 2 + /3 2 )?/> 2 V>o] ,R) 

- ( [A 2 tpQ + (a + i(3) 2 ^o + 2A(a + , R) 

- [A(i,lR 2 ) + (a-iP)(iP i' U R 2 )] 

- [2A(V 2 , \R\ 2 ) + 2(a + i/3)(VoVi, l^| 2 >] + (V>o, |i?| 2 i?)) , 



Error e (i?,i?,a) = -A -1 Rc (-([2A 2 ?/>o + 4A(a)^ 2 -0i + 2(a 2 + f3 2 )iplip ] , R) 

- < [A 2 ^l + (a + i/3) Vi^d + 2A(a + , 
[A(^ 2 ,i? 2 ) + («- i /3)(VoV'i,i? 2 )] 

- [2A(^ 2 , |i?| 2 > + 2(a + i/3)<^Vi, l^l 2 )] + (V'o, \R\ 2 R)) ■ 
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Appendix B. Dispersive Estimates 

In this section, we follow closely the work |39j on wave operators for Schrodigner operators defined on 
K. For proofs and further exposition see [39] and the references contained within. First of all, let us define 
Hq = —A and H = —A + V with the constraints on V to be discussed in the sequel. 

The wave operators, W± arc defined by 

(B.l) W ± = lim e ltH e- itH °. 

t — >oo 

Similarly, their adjoints arc defined by 

(B.2) W± = lim e itHo e- ltH P c , 

t — >oo 

where P c is the projection onto the continuous spectrum of H. The notion of wave operators is intimately 
related to the idea of distorted Fourier bases, which are discussed in detail in [5] , [TB] , [12] ■ In one dimension, 
this is directly related to the Jost solutions. These objects are studied in general in [29] and generalized 
to even a certain class of non-self-adjoint operators in [24] , 

We define a space the space of all complex-valued measurable functions <f> defined on R. such that 
(B.3) U\\ L i = J \cf>(x)\(l + \x\ydx< oo. 

Also, take the space W k ' p to be the standard Sobolev space defined by having k derivatives bounded in 
the LP norm. Then, we have the following 

Theorem 2 (Weder). Suppose that V <G L 1 for 7 > § and that for some k = 1,2, ... , V"' G L 1 /< 



or 



I = 0, 1, 2, . . . , k — 1. Then W± and W± originally defined on Wk, p fl L 2 , 1 < p < oo, have extensions to 
bounded operators on Wk,p, 1 < p < oo. Moreover, there are constants C p such that: 

(B.4) ||W±/|| W *,, < C p \\f\\ wk , P , \\W* ± f\\ wk , P < C p \\f\\ W H, P , f € W k ' p n L 2 , 1< p < oo. 

Note, there are specific requirements on the potential V which allow Theorem [2] to be extended to the 
cases p = 1 and p = oo, however we will not discuss them here. 

An important property of wave operators is that for any Borel function /, we have 

(B.5) f(H)P c = W±f(H )Wi, f(H ) = W* ± f{H)P c W ± . 

Hence, we have 

(B.6) \\e lHt P c f\\ L v = \\W ± e ltHa W* ± f\\ LP 

and using standard dispersive estimates for the linear Schrodinger operator (see for instance |35j for a 
concise overview) arrive at 

(B.7) \\e iHt P c f\\ L v < C p t-^-^\\f\\ wKP . 

Define a Strichartz pair (q, r) to be admissible if 

(B.8) 211 
K ' q 2 r 

with 2 < r < oo. Then, we arrive at the celebrated Strichartz estimates 
(B.9) \\e tHt P c uo\\ L ,w* r < \\u \\ wk ,z 

and, using duality techniques and once again the boundcdness of the wave operators, we have 

(B.IO) || f e m{t - s) P c fds\\ LqWk , r < \\f{x, t)\\ Lfw ^>, 
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where (q,r) and (q, r) satisfy (|B.8[) . . 

As a side note, using positive commutators and well crafted local smoothing spaces, from |25j . we have 
the full Strichartz estimate 

(B.H) || [* e iH ^P c fds\\ L oo L 2<\\f(x,t)\\ Lf , LS! , 

Jo * x 

where (p, q) is any allowable pair as in (|B.8[) . Now, implementing the boundedncss of wave operators on 

W k ' p spaces from [35], we have the following useful relation 

(B.12) || f e iH ^P c fds\\ Lxm < \\f(x,t)\\ L P' w ^>, 

Jo * 

where again (ft, q) is a Strichartz pair as in (|B.8[) without first going through the dispersive estimates. 

Note, as mentioned in the introduction the discussion above may be extended to the case of V having 
delta function type singularities using formalism discussed in [5]. 
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